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Abstract 

In  this  paper  we  present  a  method  for  the  design  of  analog-to- 
digital  converters  (ADCs).  This  method  computes  the  sizes  of 
the  different  components  (transistors,  capacitors,  etc.)  in  a  pre¬ 
defined  ADC  topology  so  that  the  design  specifications  are  met 
in  the  desired  process  technology. 

The  method  is  based  on  formulating  the  ADC  design  con¬ 
straints  such  as  specifications  on  power,  signal-to-noise  ratio 
(SNR),  area,  and  sampling  frequency  in  special  convex  form  in 
terms  of  the  component  sizes  of  the  ADC  and  intermediate  de¬ 
sign  variables.  More  specifically,  we  cast  the  problem  of  sizing 
the  components  of  the  ADC  as  a  geometric  program.  There¬ 
fore,  all  design  constraints  are  formulated  as  posynomial  in¬ 
equality  or  monomial  equality  constraints.  Very  efficient  nu¬ 
merical  algorithms  are  then  used  to  solve  the  resulting  geomet¬ 
ric  program  and  to  compute  the  component  sizes  of  an  ADC 
that  meets  the  desired  specifications.  The  synthesis  method  is 
fast,  and  determines  the  globally  optimal  design;  in  particular 
the  final  solution  is  completely  independent  of  the  starting  point 
(which  can  even  be  infeasible),  and  infeasible  specifications  are 
unambiguously  detected. 

This  paper  introduces  the  concept  of  hierarchical  problem 
formulation  within  a  geometric  programming  framework.  This 
modular  formulation  allows  a  high  re-use  of  the  ADC  posyno¬ 
mial  model. 

1  Introduction 

Over  the  last  five  years,  we  have  seen  a  consistent  growth  in  the 
mixed-signal  system-on-chip  (SOC)  market.  Technical  advances  in 
integrated  circuit  (IC)  manufacturing  processes  have  made  it  possi¬ 
ble  for  true  electronic  systems  [1],  such  as  cameras  and  radio  sys¬ 
tems,  to  be  integrated  in  a  single  silicon  substrate.  Since  these  elec¬ 
tronic  systems  need  an  interface  between  their  digital  components 
and  the  real  world,  there  is  a  need  for  analog  interface  circuits  to  be 
integrated  on  the  same  die  with  the  digital  components. 

The  fact  that  analog  and  digital  circuitry  have  to  co-exist  in  a 
single  substrate  has  effectively  shortened  the  required  design  time 
for  the  analog  circuitry.  When  a  new  process  technology  is  avail¬ 
able,  digital  circuitry  can  be  ported  to  it  quickly  with  the  help  of  so¬ 
phisticated  computer-aided  design  tools.  However,  since  the  analog 
part  is  designed  manually,  it  is  ported  very  slowly.  In  fact,  a  simple 
technology  port  for  an  analog  circuit  can  mean  a  complete  redesign. 

In  the  past  twenty  years,  there  has  been  extensive  research  on 


the  area  of  analog  design  automation  (see  [2]).  Approaches  to  ana¬ 
log  design  automation  can  be  classified  in  three  groups. 

•  Simulation  based  methods.  These  methods  evaluate  the  per¬ 
formance  of  the  circuit  with  a  circuit  simulator  like  SPICE 
and  search  the  design  space  using  different  types  of  optimiza¬ 
tion  algorithms.  For  example,  some  methods  use  simulated 
annealing  ASTRX/OBLX  [3],  others  gradient  search  DE- 
LIGHT.SPICE  [4],  others  a  combination  of  different  search 
methods  (MAELSTROM  [5]),  i.e.The  main  drawback  of  this 
approach  is  the  long  time  associated  with  some  simulations. 

•  Knowledge  based  methods.  This  methods  encapsulate  the 
designers  knowledge  in  some  form  of  design  plan.  Some  of 
the  most  widely  known  are  based  on  using  special  heuristics 
(like  IDAC  [6,  7]  and  OASYS  [8]  but  other  expert  systems 
have  also  been  used  (like  genetic  algorithms  SEAS  [9]).  The 
main  disadvantage  of  this  method  is  the  long  time  needed  to 
set-up  the  correct  heuristic. 

•  Equation  based  methods.  In  these  methods  the  circuit  per¬ 
formance  is  described  with  some  sort  of  analytical  design 
equations.  The  circuit  problem  is  then  cast  as  an  optimiza¬ 
tion  problem,  which  is  then  solved  by  a  numerical  algorithm. 
Some  implementation  examples  include  OPASYN  [10],  OP- 
TIMAN  [1 1],  The  largest  drawback  of  these  methods  is  their 
inaccuracy  (since  the  analytical  models  tend  to  be  too  sim¬ 
ple).  A  special  case  of  equation  based  methods  is  a  method 
based  on  formulating  the  problem  as  a  geometric  program. 
In  GPCAD  [12]  and  [13],  the  authors  describe  how  CMOS 
op-amps  and  RF  circuits  can  be  modeled  in  posynomial  form 
with  a  high  level  of  accuracy.  Since  a  geometric  program  can 
be  cast  as  a  convex  optimization  problem,  it  can  be  solved 
globally  in  a  a  very  short  time. 

The  vast  majority  of  previous  analog  design  automation  meth¬ 
ods  have  been  tested  in  small  size  circuits,  such  as  op-amps.  How¬ 
ever,  in  order  to  meet  SOC  industry  demands,  the  design  of  mid 
size  blocks  such  as  data  converters  need  to  be  automated. 

Here  we  present  a  method  for  the  design  of  a  pipeline  ADC. 
A  pipeline  ADC  has  many  more  design  variables  than  an  op-amp 
(several  hundred  versus  several  tens).  Although,  there  has  been 
some  work  on  automating  filter  design  or  ADC  design  [14],  the  dif¬ 
ference  here  is  that  we  will  size  all  design  variables  simultaneously 
by  solving  a  single  geometric  program.  In  other  works,  we  don’t 
make  decisions  at  the  system  level  and  then  design  the  circuits  at 
the  transistor  level;  we  simultaneously  decide  on  system  level  and 
transistor  level  variables.  Even  though,  we  solve  just  one  large  geo¬ 
metric  program,  the  ADC  design  is  posed  in  a  hierarchical  manner 
that  allows  re-use  of  the  formulation  when  different  building  blocks 
are  used. 

The  contribution  of  this  paper  is  to  show  that  the  design  of  a 
mixed-signal  system  composed  of  several  building  blocks  can  be 
formulated  in  an  efficient  manner  as  a  geometric  program.  The 
modular  formulation  presented  here  allows  to  effectively  model  a 
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mixed-signal  circuit.  The  fact  that  we  can  then  solve  the  problem 
globally  (rather  than  first  at  the  system  level  and  then  at  the  transis¬ 
tor  level)  results  in  much  more  optimal  designs. 

The  paper  is  organized  as  follows.  In  §2,  we  describe  the  cur¬ 
rent  custom  design  methodology  for  ADCs.  In  §3,  we  describe 
geometric  programming,  the  optimization  problem  which  is  the  ba¬ 
sis  of  the  method.  In  §4,  we  describe  the  geometric  programming 
hierarchical  formulation  of  the  ADC  design  problem.  In  §5,  we  de¬ 
scribe  how  to  use  the  method  to  design  a  specific  pipeline  ADC.  We 
start  by  introducing  the  ADC  architecture  used,  then  we  describe 
the  choice  of  design  variables  and  finally  we  present  the  posyno- 
mial  models  for  the  performance  metrics.  In  §6,  we  give  design 
examples  for  the  different  ADCs  and  show  some  design  trade-off 
curves.  In  §7,  we  end  up  with  some  conclusions  and  ideas  on  how 
to  extend  the  method. 

2  Traditional  custom  ADC  design 

In  a  traditional  custom  design  flow,  a  designer  begins  with  the  spec¬ 
ifications  for  the  ADC  that  he  needs  to  to  design.  He  starts  by 
choosing  a  suitable  architecture  or  topology.  After  that,  the  next 
step  is  component  sizing,  in  which  the  designer  determines  the  sizes 
or  values  of  the  components  for  a  given  topology  or  architecture 
that  achieve  the  requirements  or  specifications  on  the  performance 
indices.  Even  though  the  numbers  of  design  variables  and  perfor¬ 
mance  constraints  is  often  ’’small”  by  digital  circuit  design  stan¬ 
dards  (a  few  hundreds),  this  task  can  be  very  challenging,  since  in 
most  cases  all  of  the  performance  indices  are  affected  by  all  of  the 
design  variables. 

Component  sizing  in  an  ADC  is  typically  done  in  the  following 
manner.  First,  design  choices  are  made  at  the  system  level  or  top 
level.  For  example,  the  designer  first  chooses  how  many  stages  and 
how  many  bits  per  stage  a  pipeline  ADC  should  have  and  how  to 
distribute  the  overall  power  budget  amongst  stages.  Once  design 
choices  are  made  at  the  system  level,  the  designer  drills  one  level 
down  and  proceeds  to  make  additional  design  choices.  For  exam¬ 
ple,  in  the  pipeline  ADC  example,  he  will  choose  what  type  of 
amplifier  and  comparator  he  is  going  to  use  and  he  will  decide  how 
to  split  the  power  budget  of  a  stage  within  the  stage.  The  designer 
keeps  drilling  down  until  he  has  to  choose  component  dimensions 
and  values  of  bias  voltages  and  currents. 

At  each  level  of  hierarchy,  the  design  choices  are  made  in  the 
following  manner.  In  general,  the  first  step  is  to  create  a  simple 
mathematical  model  for  the  circuit.  This  model  can  be  written  in  a 
language  such  as  MATLAB  [15]  or  it  can  simply  be  a  hand  model 
and  it  is  used  to  provide  a  starting  design  point.  The  second  step  is 
to  verify  the  initial  design  point  with  a  simulator.  At  the  transistor 
level  this  is  typically  done  using  a  simulator  tool  such  as  SPICE. 
Typically,  each  building  block  is  designed  separately  (by  running  a 
large  number  of  simulations)  and  the  overall  system  is  maybe  ( only 
maybe)  simulated  once  because  of  the  very  long  simulation  times. 
Unfortunately,  this  design  methodology  is  not  optimal. 

It  is  important  to  notice  that  the  designer  does  not  simultane¬ 
ously  design  all  stages  but  rather  makes  a  lot  of  choices.  For  ex¬ 
ample,  if  he  has  a  lOOmW  power  budget  for  a  10  bit  pipeline  ADC 
he  may  decide  to  spend  20mW  in  the  first  stage,  1 5mW  in  the  sec¬ 
ond  stage  and  8.1mW  in  the  remaining  stages.  This  choice  may 
limit  the  performance  he  gets.  It  could  be  that  if  he  had  taken  into 
account  the  limitations  of  the  circuits  he  can  build,  he  would  have 
selected  30mW  in  the  first  stage,  22mW  in  the  second  stage  and 
6mW  thereafter.  It  is  hard  to  know  how  to  split  the  power  without 
knowing  the  power/performance  tradeoff  in  the  op-amps. 


3  Geometric  Programming 

Geometric  programming  (GP)  is  a  special  type  of  convex  optimiza¬ 
tion  problem  (see  [16]).  To  describe  geometric  programs,  we  first 
introduce  two  functions. 

Let  a:  be  a  vector  (*i, . . . ,  xn)  of  n  real,  positive  variables.  A 
function  /  is  called  a  posynomial  function  of  x  if  it  has  the  form 

t 

f(xi,...,xn)  =  ^CkX°lkX22k 

k=  1 

where  Cj  >  0  and  a,j  e  R.  When  there  is  only  one  term  in  the 
sum,  i.e.,  t  =  1,  we  call  /  a  monomial  function.  Note  that  posyn- 
omials  are  closed  under  addition,  multiplication,  and  nonnegative 
scaling.  Monomials  are  closed  under  multiplication  and  division. 

A  geometric  program  is  an  optimization  problem  of  the  form 

minimize  fo(x) 
subject  to  fi(x )  <1,  i  =  1, . . . 

9i(x)  =  1,  i  =  1,... 

Xi  >  0,  i  =  1, . . . ,  n, 

where  /o , . . . ,  fm  are  posynomial  functions  and  <?i , . . . ,  gp  are  mono¬ 
mial  functions. 

If  /  is  a  posynomial  and  g  is  a  monomial,  then  the  constraint 
f(x)  <  g(x)  can  be  handled  by  expressing  it  as  f(x)/g(x)  <  1. 

In  a  similar  way  if  gi  and  92  are  both  monomial  functions,  then  we 
can  handle  the  equality  constraint  gi(x)  =  gi{x)  by  expressing  it 
as 9t(x)/g2{x )  =  1. 

A  geometric  program  can  be  reformulated  as  a  convex  opti¬ 
mization  problem,  by  changing  variables  (y,  =  log  xi)  and  consid¬ 
ering  the  logs  of  the  functions  involved. 

There  are  several  methods  for  solving  geometric  programs.  One 
option  is  to  solve  the  exponential  form  of  the  geometric  program 
using  a  general  purpose  optimization  code  such  as  NPSOL  or  MI¬ 
NOS.  These  general  purpose  codes  will  in  principle  find  the  glob¬ 
ally  optimal  solution,  but  codes  specifically  designed  for  solving 
geometric  programs  offer  greater  computational  efficiency  [18].  Re¬ 
cently,  Kortanek  et  al.  have  shown  how  the  most  sophisticated  primal- 
dual  interior-point  methods  used  in  linear  programming  can  be  ex¬ 
tended  to  GP,  resulting  in  an  algorithm  with  efficiency  approaching 
that  of  current  interior-point  linear  programming  solvers  [19].  We 
use  a  simple  primal  barrier  method,  which  is  described  in  [17], 

4  Design  methodology  for  ADC  modeling 

In  this  section  we  describe  the  methodology  for  modeling  the  de¬ 
sign  of  a  pipeline  analog-to-digital  converter  as  a  geometric  pro¬ 
gram. 

The  design  constraints  are  formulated  hierarchically.  First,  the 
system  level  design  constraints  of  the  ADC  are  formulated  in  terms 
of  the  input  specifications  of  the  stages  and  design  variables  at  the 
system  level  (see  §5.2  for  more  detail).  For  example,  the  total  SNR 
is  written  as  a  posynomial  function  of  the  input-referred  noise  and 
gain  of  each  stage.  Second,  the  design  constraints  of  the  stages 
of  the  ADC  are  formulated  as  a  function  of  their  input  and  output 
specifications  and  their  design  variables.  For  example,  the  input- 
referred  noise  of  a  stage  is  written  as  a  posynomial  function  of  the 
input-referred  noise  of  the  op-amp  and  kT/C  noise  of  the  stage 
capacitors.  Finally,  the  design  constraints  of  all  stage  components 
(op-amp,  comparator,  i.e.)  are  written  as  posynomial  functions  of 
their  design  variables.  For  example,  the  input-referred  noise  of  an 
op-amp  is  described  with  a  posynomial  function  in  terms  of  tran¬ 
sistor  dimensions. 

This  hierarchical  formulation  results  in  a  modular  description 
for  the  geometric  program.  It  results  in  better  maintainability  of 
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the  implementation  and  enables  re-use  of  code  when  implementing 
the  method  for  different  ADC  topologies.  For  example,  if  we  want 
to  use  a  different  op-amp  in  the  ADC  only  the  much  smaller  op- 
amp  module  of  the  code  needs  to  be  updated. 

In  summary,  the  methodology  consists  of  the  following  steps: 

1.  Selecting  ADC  topology.  Depending  on  the  application  of  the 
ADC  a  certain  topology  for  the  ADC  is  chosen. 

2.  Defining  levels  of  hierarchy  and  corresponding  input ,  output 
and  design  variables.  As  we  will  describe  in  §  5.2,  the  ADC 
is  divided  into  three  levels  of  hierarchy.  At  each  level  we 
define: 

•  Input  variables.  These  are  the  input  specifications  to  a 
circuit  block.  For  example,  in  the  top  level  the  input 
variables  are  just  the  input  specifications  to  the  con¬ 
verter  like  the  SNR  specification. 

•  Design  variables.  These  are  the  variables  that  are  com¬ 
puted  at  that  level  of  hierarchy.  For  example,  in  the  top 
level  a  design  variable  would  be  the  number  of  stages 
the  ADC  needs. 

•  Output  variables.  These  are  the  specifications  imposed 
in  the  circuits  a  level  below.  For  example,  in  the  top 
level  the  required  noise  level  of  each  stage  would  be  an 
output  variable. 

Each  hierarchy  level  has  a  minimal  number  of  defining  input, 
output  and  design  variables  that  are  sufficient  to  describe  the 
behavior  of  the  ADC  and  the  interaction  between  levels. 


per  stage  resolution  stage  is  as  follows.  The  analog  input  is  sam¬ 
pled  into  capacitors  Ci  and  Ci  during  the  sampling  phase.  During 
the  transfer  phase,  the  ADC  performs  a  coarse  quantization  of  the 
input  signal  which  is  subtracted  from  the  held  signal  and  then  am¬ 
plified.  The  residue  of  the  subtraction/amplification  operation  is 
passed  down  to  the  next  stage  for  fine  quantization.  Much  more  de¬ 
tail  about  the  basic  operation  of  pipeline  ADC  can  be  found  in  [20]. 


Figure  1 :  Single-ended  implementation  of  a  pipeline  ADC  stage 


Although  we  have  limited  ourselves  to  same  number  of  bits  per 
stage,  the  method  can  be  extended  to  architectures  composed  of 
stages  with  different  bits  of  resolution  (e.g.,  first  stages  two  bits 
and  last  stages  one  bit). 


3.  Writing  (posynomial)  ADC  design  equations  at  each  hierar¬ 
chy  level  in  terms  of  input,  output  and  design  variables  of 
sub-blocks.  At  this  step,  the  design  constraints  of  the  ADC 
are  put  in  posynomial  inequality  form  in  terms  of  the  input, 
output  and  design  variables  of  the  sub-blocks.  This  step  and 
the  previous  step  introduce  a  hierarchical  methodology  for 
writing  ADC  design  equations  in  terms  of  the  ADC  compo¬ 
nent  sizes. 

4.  Formulate  problem  as  geometric  program  and  solve  for  com¬ 
ponent  sizes  using  numerical  algorithm.  Once  design  con¬ 
straints  are  put  in  posynomial  form,  the  ADC  design  problem 
is  cast  as  a  geometric  program  and  hence  it  can  be  readily 
solved  using  efficient  numerical  algorithms. 

Note  that  even  though  the  formulation  is  hierarchical,  the  re¬ 
sulting  geometric  program  is  solved  in  a  flat  manner.  In  other 
words,  design  equations  at  all  levels  are  solved  simultane¬ 
ously.  This  can  result  in  a  very  large  optimization  problem 
but  given  the  efficiency  of  geometric  program  solvers  it  does 
not  become  an  issue. 


5.2  Design  and  I/O  variables 

In  this  section  we  show  how  a  variety  of  performance  measures  and 
constraints  can  be  formulated  using  geometric  programming.  We 
define  three  hierarchy  levels  (see  Figure  2), 


Input  specs 


Top  level 


-  Number  of  stages 


1  Stage  specs 


Stage  level 


-  Capacitors 


Circuit  specs 


Transistor  level 


-Transistor  sizes,  bias  values 


5  Posynomial  model  of  a  pipeline  ADC 
5.1  Pipeline  ADC  architecture 

To  simplify  the  discussion  we  consider  a  specific  op-amp  topology, 
the  1  bit  per  stage  ADC  of  Figure  1.  In  practice  the  1  bit  per  stage 
ADC  is  implemented  with  a  1.5  bit  per  stage  but  the  assumption  of 
1  bit  per  stage  simplifies  the  explanation  of  the  method. 

Pipeline  ADCs  consist  of  a  set  of  stages  connected  in  series. 
Figure  1  shows  a  single-ended  implementation  of  a  one  bit  per  stage 
resolution  pipeline  ADC  stage.  Each  stage  consists  of  four  capac¬ 
itors  (two  in  a  single  ended  implementation),  a  digital  to  analog 
converter  (DAC),  an  operational  amplifier  and  ten  switches  (six  in 
a  single  ended  implementation).  The  basic  operation  of  the  one  bit 


Figure  2:  Hierarchical  design  of  a  pipeline  ADC 

•  Hierarchy  level  one:  System  level.  The  input  specifications 
to  this  level  are  the  input  specifications  of  the  converter: 

-  Bits  of  resolution 

-  Power  consumption 

-  Input  signal  bandwidth 

-  Maximum  sampling  frequency 

-  Clock  characteristics 

-  Supply  voltage 
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*  Dimensions  of  bounding  area. 


-  Input  range  (Kef) 

-  Input  signal  common  mode 

-  Dynamic  range  (DR) 

-  Linearities  (INL  and  DNL ) 

-  Signal  to  noise  plus  distortion  ratio  ( SNDR ). 

-  Dimensions  of  bounding  area 

The  design  variables  are 

-  Number  of  stages 
The  output  specifications  are 


•  Hierarchy  level  three:  Circuit  level.  The  circuit  level  is  the 
third  and  last  level  of  hierarchy.  The  input  specifications  are 
the  same  as  the  output  specifications  of  the  previous  level. 

The  design  variables  are 

-  Width  of  each  transistor 

-  Length  of  each  transistor 

-  Number  of  fingers  of  each  transistor 

-  Values  of  passive  components  (capacitors  and  resistors) 

-  Value  of  bias  currents  and  bias  voltages. 


-  Single  pipeline  stage  noise 

-  Single  pipeline  stage  power 

-  Single  pipeline  stage  area 

-  Single  pipeline  stage  offset 

-  Single  pipeline  stage  nonlinearities 

-  Specifications  for  clock  generation  circuitry,  bias  cir¬ 
cuitry  and  voltage  reference  generation  circuitry. 

-  Dimensions  of  bounding  area  for  each  stage 

We  focus  our  attention  on  the  pipelined  stages.  The  addi¬ 
tional  circuitry  (clock  generation,  bias,  and  voltage  reference 
generation)  will  be  handled  with  detail  in  the  circuit  imple¬ 
mentation  phase  of  this  project. 

•  Hierarchy  level  two:  Stage  level.  The  second  level  of  hier¬ 
archy  is  composed  of  the  building  stages.  In  this  level,  the 
input  specifications  are  the  same  as  the  output  specifications 
of  the  top  level.  The  design  variables  are  the  capacitor  sizes. 

The  output  specifications  are 

-  Op-amp  specifications, 

*  Unity-gain  bandwidth 

*  Settling  time 

*  Load  capacitance 

*  Power 

*  Gain 

*  Slew  rate 

*  Output  swing  and  common  mode  level 

*  Area 

*  Nonlinearities 

*  Input-referred  noise 

*  Offset  voltage 

*  Dimensions  of  bounding  area 

-  Switches  specifications 

*  On  resistance 

*  Turn-on  time 

*  Area 

*  Nonlinearities 

*  Dimensions  of  bounding  area 

-  Comparator  specifications 

*  Power 

*  Area 

*  Speed 

*  Offset  voltage  and  meta-stability  range 

*  Input-referred  noise 


5.3  Performance  metrics 

1.  System  level  design  constraints 


•  If  we  assume  equal  resolution  per  stage,  the  number  of 
stages  M  is  defined  by  the  posynomial, 


M  =  N/B,  (2) 

where  N  is  the  number  of  bits  of  resolution  of  the  con¬ 
verter  and  B  is  the  number  of  bits  of  resolution  per 
stage.  Note  that  equation  2  is  a  monomial. 

In  the  development  of  this  CAD  tool,  we  have  assumed 
that  the  bits  per  resolution  per  stage  (B)  are  equal  in 
all  stages.  B  is  an  input  to  the  tool  and  not  a  design 
variable.  Therefore,  in  order  to  evaluate  whether  one, 
two  or  three  bits  of  resolution  per  stage  is  better,  one 
must  run  the  CAD  tool  three  times  with  three  different 
topology  specifications  (B  =  1,  B  =  2  and  B  =  3). 

•  The  dynamic  range  is  given  by  (see  [20]), 


DR  =  10  log 


(2JV~1  A)2/2 


(3) 


where  np  ,  the  noise  power  at  the  input  of  the  converter, 
is  given  by 


M 

np  =  nQ  + 

i=l 


Tlstagei 

GW”1)  ’ 


(4) 


where  nstagei  is  the  input  referred  noise  of  the  j'th  stage 
and  riQ  is  the  quantization  noise  given  by 


_  (Vret/2N)2 

12  12 


(5) 


Therefore  if  we  want  to  impose  a  condition  on  a  max¬ 
imum  allowed  dynamic  range,  we  would  impose  the 
following  posynomial  constraint, 


M 


Tl'tase;  „  (2JV~1A)2/2 


G20- 1) 


10- 


(6) 


The  SNDR  is  given  by 


SNDR  =  10  log 


(2jv"1A)2/2 

Tip  — b  hp 


(7) 


where  hp,  the  harmonic  power  at  the  input  of  the  con¬ 
verter,  is  given  by. 


4 


(8) 


••'stagei 

GW-1) ' 


Therefore  if  we  want  to  impose  a  condition  on  a  max¬ 
imum  allowed  SNDR,  we  would  impose  the  following 
posynomial  constraint, 


M  M 

"«+£SS+E 


8  =  1 


•  f'Stagei 

GW-1) 


(2jv-1A)2/2 


10 — ro 


(9) 

•  The  power  of  the  converter  is  given  by  the  posynomial, 


p=j2Pi+p„ 


(10) 


where  Pi  is  the  power  of  the  ith  stage  and  Pv  is  the 
power  consumed  by  the  clock  generation  circuitry,  bias 
circuitry,  voltage  reference  generation  circuitry  and  buffers. 

•  The  sampling  frequency  translates  into  the  following 
monomial  constraint, 

f ‘sampling  ^  f stage, i-  (11) 

This  condition  imposes  a  minimum  operating  frequency 
for  each  stage. 

•  The  area  of  the  converter  is  given  by  the  posynomial, 

M 

A  —  ^  ^  Ai  +  Ajj  +  -Aroutej  (12) 

i=l 

where  Ai  is  the  area  of  the  ith  stage,  Ar/  is  the  area 
of  the  clock  generation  circuitry,  bias  circuitry,  voltage 
reference  generation  circuitry,  buffers  and  decoupling 
capacitors  and  Aroutc  is  the  area  due  to  the  routing  of 
each  stage  and  additional  circuitry. 

•  The  floorplan  of  the  converter  can  also  be  taken  into 
account.  For  example,  consider  the  simple  transistor 
layout  of  Figure  3,  where  a  bounding  box  defined  by 
dimensions  Xd  and  yd  is  given  to  lay  out  a  10  bit  ADC. 
This  ADC  is  composed  of  ten  stages,  digital  circuitry 
and  bias  circuitry.  Assuming  we  have  the  fixed  floor- 
plan  of  Figure  3  we  add  constraints  to  make  sure  that 
all  building  blocks  fit  in  the  bounding  box.  For  this 
simple  layout,  the  constraints  would  be 
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Figure  3:  Floorplan  of  a  10  bit  ADC 


•  The  output  of  the  gain  block  is  given  by 


V(n  -  1)  = 


[(‘  +  §)  V(n)  +  i  (-D<">vtcflP  -  , 


03) 


where  Ao  is  the  op-amp  open-loop  gain  and  r  is  the 
closed-loop  time  constant.  This  time  constant  is  related 
to  the  unity-gain  bandwidth  ( UGBW )  and  open-loop 
pole  (pi)  of  the  op-amp  by. 


r= _ tlA _ ~ _ 1 _ (i+CA 

pi  (1  +  £  +  Ao)  2tt UGBW  \  +  C2)’ 

(14) 

In  order  to  operate  at  the  required  frequency,  the  fol¬ 
lowing  posynomial  constraint  must  be  imposed  for  the 
first  stage, 


Tsettle  +  Tslew  +  Tclk  < 


2/samplii 


(15) 


where  Tscttie  =  2 t(N  +1),  Tcu,  is  the  delay  due  to  the 
clocking  scheme  and  the  worst  value  for  Tsicw  is  given 

by. 


Tslew 


AVin 


SR 


Vref 

2  SR' 


(16) 


Note  that  N  is  full  precision  for  the  first  stage,  N/Gi 
for  the  second  stage,  N/(GiG2)  for  the  third  stage  and 
so  on.  This  means  that  the  settling  time  constraint  (15 
is  more  constraining  for  the  first  stages. 

Note  also  that  Tsettie  is  a  posynomial,  Tc ik  is  an  input 
specification  and  therefore  a  constant  and  that  Tsiew  is 
a  posynomial.  This  means  that  equation  (15)  is  also  a 
posynomial. 

•  The  load  capacitance  of  the  operational  amplifier  in  stage 
ith  (CL;)  during  the  charge  transfer  phase  is  given  by, 


2.  Stage  level  design  constraints 


CLi  —  Ctransfer)Stagejj  -f-  CsamplingjStage]j_i  “b  Gp  — 

(Cjnput  “b  Cl ,j)  C2,i 


Cjnput  +  CU  +  C2,i 


+  {C2,  i+1  +  Cl.i+i  +  C'comp)  *4*  Cp , 

(17) 
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where  Cinput  is  the  input  capacitance  of  the  op-amp 
in  stage  ith,  Cc omP  is  the  input  capacitance  of  the  com¬ 
parator  in  the  i+lth  stage  and  Cp  is  the  parasitic  capac¬ 
itance  from  the  switches.Note  that  equation  (17)  is  not 
posynomial.  Even  though  this  equation  is  not  posyno- 
mial  it  can  be  very  well  approximated  when  the  bits  per 
resolution  per  stage  is  known.  For  example  if  B  =  1, 
Ci  =  C2  and  we  can  write  the  posynomial, 


+  (<72  ,*+l  +  Ci  ,t'+l  4”  ^comp  )+cm) 

•  The  power  of  the  stage  is  given  by  the  posynomial 

jP  —  Pop— amp  +  ■f^omparator  j  (19) 

where  Pop-amp  is  the  power  consumed  by  the  op-amp 
and  ^comparator  is  the  power  consumed  by  the  com¬ 
parator. 

•  The  stage  gain  is  given  by 

Gain  =  1  +  ^-  =  2b  (20) 

G2 

where  B  is  the  bits  per  resolution  per  stage.  Since  we 
fix  B  at  the  beginning  of  the  problem,  we  can  write  the 
following  monomial  constraint, 

=  2fl  -  1  (21) 

C/2 

•  The  area  of  a  stage  is  given  by  the  posynomial  equation, 


A  —  -Aop-amp  4“  ^comparator  4"  2 Aci  4* 

2Ac2  4“  E  ■^switch, j  4"  -Aroute?  (22) 
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where  Aop-amP  is  the  op-amp  area,  ^comparator  is  the 
comparator  area,  Ac\  is  the  area  of  capacitor  Ci ,  Ac2 
is  the  area  of  capacitor  C2,  ^switch, j  is  the  area  of  each 
switch,  and  +routc  is  the  routing  area. 

•  The  noise  for  the  ith  stage  is  given  by  the  sum  of  the 
thermal  noise  in  the  switches  and  the  amplifier  noise 
(see  [21]  for  a  detailed  derivation), 


constraints  on  the  op-amp  are  given  by  the  following 
monomial  constraints, 

Tout, max  >  (1  +  a)  (VcM  + 

Vout.min  <  (1  —  a)  (VcM  —  ■^t)  , 

where  a  is  a  factor  that  accounts  for  the  linearity  mar¬ 
gin  (typically  a  few  hundred  millivolts  are  reasonable), 
Vcm  is  the  common  mode  output  voltage  and  Vref  is 
the  reference  voltage. 

•  In  order  to  minimize  the  nonlinearities  due  to  the  op- 
amp  gain  nonlinearity  and  reduce  the  finite  op-amp  gain 
effect,  the  gain  of  the  amplifier  has  to  be  sufficiently 
high.  We  can  achieve  this  by  imposing  the  monomial 
constraint, 

Ao.i  >  2Ni+1,  (26) 

where  A0,i  is  the  op-amp  gain  of  the  ith  stage  and  N, 
is  the  precision  bits  required  in  the  ith  stage. 

•  In  order  to  achieve  a  good  settling  behavior,  the  phase 
margin  for  the  amplifier  has  to  be  sufficiently  high.  In 
practice,  guaranteeing  that  the  phase  margin  is  at  least 
seventy  degrees  ensures  a  good  settling  behavior. 

•  The  causes  of  nonlinearities  in  a  stage  are  several  [22]. 
The  cause  of  DC-nonlinearity  include:  amplifier  and 
comparator  offsets,  capacitor  mismatch  and  nonlinear¬ 
ity,  amplifier  finite  gain  and  nonlinearity,  and  switch 
charge  injection.  Since  we  are  imposing  conditions  (26 
and  25),  we  can  neglect  the  nonlinearity  due  to  the  am¬ 
plifier  gain.  Capacitor  mismatch  and  nonlinearity  data 
are  process  dependent  parameters  and  as  such  they  are 
an  input  to  the  tool.  Therefore  we  can  write  (see  [22]), 

frstage  —  Vref )  +  (offsetop_amp)  + 

(offsetcomparator)  +  (^ci)  ,  (27) 

where  ec  is  the  percent  mismatch  and  nonlinearity  of 
the  capacitors,  offsetop-amp  is  the  input  offset  of  the 
op-amp,  offsetcomparator  is  the  input  offset  voltage  of 
the  comparator  and  qc\  is  the  total  charge  injection  due 
to  the  switches  of  the  stage  (see  [23]). 

Dynamic  nonlinearity  has  not  been  modeled.  In  prac¬ 
tice,  this  is  a  difficult  performance  metric  to  model  and 
simulate. 

•  The  switches  must  have  an  enough  small  associated 
time  constant. 


e?'2Kr(Jr+s)+e:--<’  <*» 


where  eLP,i 
by 


is  the  input-referred  op-amp  noise  given 


So 

At’ 


(24) 


where  So  is  the  op-amp  input-referred  white  noise  den¬ 
sity  (given  by  a  posynomial ,  see  [  1 2] )  and  r  is  the  closed 
loop  time  constant  given  by  equation  (14). 

Since  S0  is  a  posynomial  and  r  a  monomial,  equa¬ 
tion  (23)  is  a  posynomial. 

•  The  output  swing  of  the  op-amp  has  to  be  at  least  as 
large  as  the  reference  signal.  In  order  to  guarantee  high 
linearity,  we  leave  some  extra  margin.  The  output  swing 


B,  o : :  Cs  ,v  l//stage,  (28) 

where  Csw  is  the  capacitance  seen  by  each  switch.  If 
we  impose  a  factor  of  ten  difference  between  the  switch¬ 
ing  time  constant  and  the  sampling  frequency  equa¬ 
tion  (28)  is  a  posynomial. 

•  By  fixing  a  priori  the  floorplan  of  each  stage  we  can 
impose  posynomial  constraints  on  the  floorplan  (in  a 
similar  way  to  what  we  did  in  the  top  level  of  hierar¬ 
chy). 

3.  Circuit  level  design  constraints 

The  problem  of  modeling  an  op-amp  as  a  geometric  program 
has  been  described  in  previous  papers  [12]  and  will  not  be 
described  here.  Modeling  a  comparator  and  switch  opera¬ 
tion  can  be  done  in  a  similar  manner.  Recently,  new  meth¬ 
ods  for  automatically  developing  posynomial  circuit  models 
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have  been  reported  [24].  These  new  methods  can  also  be 
used  successfully  (especially  because  of  the  small  number  of 
variables  associated  with  the  comparator  and  the  switches). 

One  could  envision  defining  a  fourth  level  of  hierarchy,  de¬ 
fined  by  the  actual  transistors.  In  other  words,  in  the  third 
level  we  would  formally  write  op- amp  specifications  in  terms 
of  transistor  parameters  and  in  the  fourth  level  we  would  re¬ 
late  transistor  parameters  to  transistor  sizes.  Transistor  be¬ 
havior  can  be  modeled  in  posynomial  manner.  A  simple 
posynomial  model  ( GP1  model)  was  described  in  [12]  but 
more  complex  posynomial  models  can  also  be  developed. 

6  Design  examples 

We  have  implemented  a  simple  CAD  tool  for  the  design  of  pipeline 
ADCs  with  one  bit  per  resolution  stages.  We  have  modeled  in 
posynomial  form  a  gain  boosted  op-amp  [12],  a  dynamic  compara¬ 
tor  and  a  set  of.simple  NMOS  switches.  The  input  to  the  CAD 
tool  are  the  converter  specifications:  number  of  bits,  maximum 
SNR,  maximum  sampling  frequency,  maximum  area  and  maximum 
power.  The  output  are  the  sizes  of  the  capacitors  and  the  transistor 
sizes  for  the  op-amp,  comparator  and  switches. 

The  tool  is  written  in  C  code.  After  reading  a  set  of  speci¬ 
fications,  a  geometric  program  consisting  of  all  the  posynomial 
and  monomial  equations  shown.  For  a  12  bit  ADC,  the  problem 
contains  2364  variables  and  it  is  solved  in  under  10  minutes  on  a 
400MHz  Pentium  PC  with  128MB  of  memory. 

In  Figure  4  we  show  the  tradeoff  curve  power  versus  sampling 
frequency  for  a  12  bit  ADC  and  a  10  bit  ADC.  We  impose  an  SNR 
of  70dB  for  the  12  bit  ADC  and  of  58dB  for  the  10  bit  ADC  (in 
both  cases  a  loss  of  4dB  due  to  circuit  noise).  This  is  done  by  re¬ 
peatedly  solving  the  ADC  design  problem  (minimize  power)  while 
varying  the  sampling  frequency  specification.  We  observe  (as  ex¬ 
pected)  that  the  higher  the  required  sampling  frequency  the  more 
power  is  needed.  What  this  curve  allows  us  to  understand  is  how 
exactly  these  specifications  trade-off.  Note  that  less  power  is  re¬ 
quired  for  a  10  bit  ADC  but  it  also  has  a  smaller  SNR.  Typical 
capacitor  mismatches  limit  the  resolution  of  high  speed  ADCs  to 
10  bits  and  calibration  techniques  are  needed  to  achieve  higher  res¬ 
olution.  Here  we  have  simulated  the  use  of  calibration  by  using  a 
very  small  capacitor  mismatch  parameter. 


Minimum  power  vrosus  sampling  frequency  tor  12  bh  and  10  bit  A2D 


Figure  4:  Minimum  power  for  12  bit  and  10  bit  ADC. 


Recently,  there  has  been  some  work  done  to  define  optimum 
scaling  in  pipeline  ADC  [25].  However  very  simple  models  were 


used  to  determine  this  optimum  scaling.  Using  the  approach  pre¬ 
sented  here  we  can  quickly  determine  what  is  the  optimum  scaling. 
Figure  5  and  figure  6  show  the  optimum  power  per  stage  and  op¬ 
timum  capacitance  scaling  for  the  case  of  a  12  bit  ADC  operating 
with  a  maximum  sampling  frequency  of  20MHz.  We  observe  that 
from  the  fifth  stage  on,  all  stages  are  the  same  size  and  scaling 
only  takes  place  in  the  initial  stages.  The  reason  is  that  after  the 
fifth  stage,  the  op-amp  behavior  is  determined  by  its  own  parasitics 
rather  than  by  the  switching  capacitances  which  are  quite  small  in 
later  stages  (small  capacitances  in  later  stages  are  possible  because 
their  kT/C  noise  is  attenuated  by  the  gain  of  the  initial  stages). 
The  fact  that  the  last  stages  are  identical,  allows  to  reduce  the  size 
of  the  problem  by  almost  half.  The  idea  is  to  assume  the  first  six 
stages  are  unique  and  the  last  six  stages  are  identical. 


Power  scaling  In  12  bit  A20, 20MHz 


Figure  5:  Power  scaling  for  12  bit  ADC,  20MHz. 


Capacitance  scaling  In  12  bit  A2D,  20MHz 
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Figure  6:  Capacitor  scaling  for  12  bit  ADC,  20MHz. 


7  Conclusions  and  extensions 

In  this  work  we  have  shown  how  to  use  geometric  programming 
to  design  a  pipeline  ADC.  Since  the  method  is  very  efficient,  we 
can  simultaneously  size  all  pipeline  stages.  Since  we  a  solving  a 
convex  problem,  the  results  obtained  are  globally  optimal. 
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An  important  feature  of  geometric  programming  based  design 
is  the  ability  to  develop  robust  designs,  i.e.,  designs  that  meet  the 
required  specifications  are  met  for  a  variety  of  different  technology 
parameter  values  and  operating  conditions.  To  do  this,  one  just 
needs  to  replicate  the  design  constraints  for  the  different  scenarios, 
which  is  practical  only  because  the  computational  effort  for  solving 
geometric  programs  grows  approximately  linearly  with  the  number 
of  constraints.  The  ability  of  creating  robust  designs  is  perhaps 
the  most  important  feature  of  geometric  programming  based  design 
since  most  of  the  time,  designers  are  more  concerned  with  having 
a  circuit  that  works  over  all  comers  than  a  globally  optimal  circuit. 

The  method  presented  does  not  create  new  ADC  topologies. 
What  it  does  is  size  a  previously  defined  ADC  architecture.  How¬ 
ever,  the  modularity  of  the  method  allows  to  effectively  create  new 
topologies.  For  example,  if  we  have  a  library  of  ten  op-amps  and 
five  comparators,  we  can  envision  several  possible  combinations 
(all  stages  have  same  op-amp  and  comparator;  first  four  stages  use 
a  certain  op-amp  and  certain  comparator  and  the  rest  of  stages  use 
a  different  type,  i.e.).  In  order  to  evaluate  all  different  combinations 
one  does  not  need  to  re-formulate  the  entire  ADC  design  problem. 
One  only  needs  to  formulate  the  module  in  question  as  a  geomet¬ 
ric  program  following  the  convention  of  required  input,  output  and 
design  variables.  The  search  of  all  possible  combinations  is  a  com¬ 
binatorial  problem.  Since  geometric  programs  can  be  solved  very 
fast,  one  can  just  search  the  space  by  solving  each  possible  ADC 
architecture.  At  some  point,  however  this  search  becomes  too  long 
(when  too  many  combinations  are  possible)  so  more  efficient  search 
schemes  would  be  needed  (this  investigation  is  not  the  purpose  of 
this  paper). 

The  approach  we  have  described  is  not  limited  to  pipeline  ADCs. 
Hierarchical  GP  formulation  can  be  used  to  describe  other  data  con¬ 
verters  and  other  mixed-signal  blocks  such  as  phase-locked-loops. 
The  methodology  is  identical  to  the  one  described  in  §4. 
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